Is HBT really puzzling? 
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Abstract. Two-particle correlations from RHIC have provided a surprising snapshot of the final 
state at RHIC. In this talk I discuss the nature of the HBT puzzle and attempt to delineate several 
\ factors which might ultimately resolve the issue. 
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Correlation, or HBT, measurements provide our best insight into the space-time 
development of relativistic heavy-ion collisions [1]. Immediately after the first two- 
^ ■ pion correlations from RHIC were presented in 2001, the term "HBT puzzle" came into 

common usage. The term was inspired by the failure of some of our most sophisticated 
dynamical models of the collision. However, there is no problem in describing the data 
with simple parameterizations of the break-up space-time profile. In this talk, I will 
review one of these simple parameterizations, the blast-wave model, to illustrate the 
unsettling aspects of the inferred parameters. I then discuss several factors which might 
^ ; ultimately lead to a more satisfying interpretation of the data. 

Numerous parameterizations of the breakup space-time geometry can be applied to fit 
correlation data from RHIC. Among these are the Buda-Lund model [2] and numerous 
variations of the blast- wave model [3, 4, 5]. Minimal blast wave models incorporate 
four common parameters: the breakup temperature T, the radius of the fireball surface 
7?, the collective velocity at the surface v± 9 and the breakup time T, which is related 
to the collective- velocity gradient along the beam direction by the relation v z = z/t by 
virtue of assuming boost invariance. Additionally, chemical potentials can be added to 
normalize the spectra. Parameters for the surface diffuseness or the temporal duration of 
the emission At are often included. The Buda-Lund model also allows for a temperature 
gradient. 

Pionic observables do not tightly constrain the parameters by themselves [6], but after 
fitting the spectra of heavier particles the following parameters emerge from blast- wave 
fits: T ~ 110 MeV, v± ~ 0.7c, R ~ 12 - 13 fm, and T ~ 9 fm/c. The duration of the 
emission appears to be sudden, less than 5 fm/c, and roughly consistent with zero. The 
quality of the fits with data are illustrated in Fig. 1. 

On inspection, there are two puzzling aspects of these numbers. First, the breakup 
density is remarkably high. If one takes the 1000 hadrons per unity rapidity observed at 
RHIC and divide it by the volume V = ttiR 2 , one obtains densities higher than nuclear 
matter density, in the range of 0.18 fm -3 . Assuming cross sections of 20 mb, estimates 
of the mean free path are in the neighborhood of 2.5 fm, which is much smaller than 
the 25 fm diameter of the fireball. The second surprising aspect of the fit is that it 
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FIGURE 1. Blast wave fits to source sizes from Lisa and Retiere [3]. 



implies the fireball, whose original diameter was 6 fm, expanded by 7 fm in 10 fm/c. 
Given the surface velocity was 0.7c, this implies the surface unphysically accelerates 
instantaneously to its final velocity. 

The density, though surprising, can be readily explained. First, due to collective 
expansion, the density will have fallen substantially by the time particles move a few 
fm. Secondly, at the time of breakup, the expansion has become isotropic. For non- 
relativistic Gaussian-like expansions of particles of the same mass, collisions do not alter 
the outgoing phase space density, and thus cease to have an effect on the overall source 
size [7]. Thus, the phase space distribution tends to effectively freeze out before particles 
have their last collisions. In fact, microscopic simulations incorporating cross sections of 
tens of mb often yield source sizes smaller than those measured at RHIC. The freezing- 
out of the source size can be thought of as a competition between the cooling, which 
pushes the system toward smaller source sizes and the growth of the overall source. For 
the qualifiers mentioned above, these two balance perfectly. In practice, collisions of 
pions with heavier particles tend to reduce the pion source and increase the apparent 
source size of heavier particles such as protons, as entropy moves from pions to protons. 
In fact, source sizes from microscopic models usually vary remarkably little as a function 
of the assumed cross section [8]. 

Source sizes are constrained by entropy, which can be calculated from final-state 
phase space densities, which can in turn be extracted from experiment given measured 
radii, 
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These expressions were applied in [9] to extract entropy and phase space density in 
central 130A GeV collisions at RHIC. Phase space densities at this and lower energies 
from [1] are displayed in Fig. 2. Note that the phase space densities tend to rise, then 
saturate near the top SPS energies. At this point the phase space density, which is 
approximately frozen after chemical freezeout in an isentropic expasion, saturates when 
the entirety of the fireball surpasses temperatures of 170 MeV [10], the temperature for 
which particle ratios suggest chemical compositions have frozen out [11]. 

Pionic source volumes are a factor of two smaller than those predicted by hydrody- 
namic models based on lattice-gauge-theory-inspired equations of state. This represents 
a deficit of ln(2) units of entropy per pion, which might suggest that the equations of 
state used for the hydrodynamic calculations in [12] is excluded. However, an analysis 
of the entropy carried by all the particles is consistent with entropy expected from lattice 
calculations [9] as shown in Fig. 2. Since baryons tend to carry twice as much entropy 
per particle as pions, a modest increase in the population of baryons is able to absorb the 
~ 20% shortfall in pionic entropy. If one knows s(e) 9 the entropy density as a function 
of the energy density, one can also derive the pressure P(e) from thermodynamic identi- 
ties. A higher-entropy equation of state, i.e., one where the effective number of degrees 
of freedom is higher, must also have lower pressure than an equation of state whose en- 
tropy rises more slowly with energy. Since entropy can increase during the expansion, 
the fact that the measured entropy is consistent with the lattice equation of state places 
a lower-bound on the stiffness of the equation of state. A very stiff equation of state, 
like the pion gas example illustrated in Fig. 2, could in principle be consistent with the 
final-state entropy, but would require the entropy to double during the expansion. Most 
mechanisms for generating entropy such as viscosity and shocks are expected to provide 
^ 10% increases in entropy. 

Unfortunately, there are no systematic studies of the dependence of source radii as 
a function of the equation of state. However, a compilation of several hydrodynamic 
and microscopic transport models presented in Fig. 3 shows a trend of stiffer equations 
of state resulting in smaller radii, as expected. By incorporating more resonances, or 
by including strings, microscopic models effectively lower the pressure. The largest 
source sizes resulted from hydrodynamic calculations using an equation of state with 
a large latent heat. Although it is premature, one might conclude that the equation of 
state appears somewhat stiffer than some of the lattice models, but not nearly so stiff 
as that of a pion gas. Taken as a whole the model calculations are encouraging as they 
illustrate the sensitivity of correlation meausurements to the equation of state. 

Even though the average source sizes can be reproduced with some of the transport 
models represented in Fig. 3, none of the models does a good job of reproducing the 
p t dependence for all three dimensions. Part of this difficulty lies in the fact that, as 
stated earlier, models tend to require some time to accelerate transversely and thus have 
difficulty reproducing the final-state geometry described by the blast-wave model. This 




FIGURE 2. Average pionic phase space densities as a function of p t (left panel) rise with increasing 
beam energy, but saturate at upper SPS energies. As a function of the initial energy density at 1 fm/c, 
expected entropies per unit rapidity can be calculated given the equation of state (right panel). A blast- 
wave model with a chemical potential of 70 MeV for pions roughly reproduces the RHIC results. The 
Bjorken estimate probably underestimates the true energy density by ~ 50%, thus the lattice equation of 
state seems consistent with the experimentally extracted dS/dy which is slightly above 4000. 



aspect of the HBT puzzle might be resolved by some combination of the following 
improvements or alterations of dynamic models: 

1. Reduced emissivity. Sequential evaporation from the surface tends to increase 
^out/^side beyond the value of unity observed. Since pions with p t > 100 MeV/c 
move faster than the surface, those pions emitted earlier get ahead of other pions 
and extend the shape of the phase space distribution. Any mechanism to reduce the 
emissivity of the surface, such as super-cooling [13], would help the system reach 
the sudden freeze-out picture of a blast-wave model. 

2. Non-infinite longitudinal extent. Many of the hydrodynamic models assume boost 
invariance. A finite extent along the beam direction lowers R\ ong , and since R\ ong is 
related to the inferred lifetime, accounting for the finiteness should lead to longer 
lifetimes which gives the matter more time to expand transversely to the large 13 
fm size. 

3. Longitudinal acceleration. Boost invariant models also neglect acceleration along 
the beam axis. Accounting for this acceleration alters the connection between the 
velocity gradient and the lifetime, dv z /dz ~ 1/t, and should result in somewhat 
longer times [14]. 

4. Shear viscosity. The high velocity gradient along the beam axis at early times can 
lead to shear effects which lower the pressure along the beam axis while raising the 
transverse pressure [15, 16]. This effect could be especially significant if the early 
stage is well described by longitudinal classical fields. For non-interacting electric 
fields along the z axis, the components of the stress energy tensor are T xx = T yy = £, 
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FIGURE 3. Dimensions of several models from identical-pion correlations compared to RHIC data as a 
function of the average pair momentum k t . A basic trend is apparent that those models with soft equations 
of state overestimate the source sizes while the stiffest example, a pion gas (MPC) grossly underestimates 
the source sizes. None of the models quantitatively describe the k t dependence of all the dimensions. 
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5. Refraction through the mean field. An attractive mean field for pions could refract 
trajectories and lead to phase space distributions whose R out dimension exceeds the 
physical size of the emitting source [17, 18, 19], especially the 7? s id e projection. 

None of these effects are expected to represent more than a 10% correction, except for 
the refractive effects at very low p t9 but they all push the analysis of the data toward 
a more physical interpretation. Thus, the HBT puzzle might well have originated from 
a conspiracy of effects, each of which borders on being negligible when considered by 
itself. 

Unraveling the space-time picture of these collisions will undoubtedly benefit from 
improved measurement, especially the measurement of other correlations besides pion 
correlations [20, 21]. However, the main impediment toward resolving correlation mea- 
surements remains the lack of a cohesive well-tested modeling framework. For instance, 
there is the need of a detailed comparison of hydrodynamic and Boltzmann approaches 
using the same effective equations of state, so that the role of viscosity, and the sensitivity 
to breakup conditions might be quantified. The random comparison of different models, 
each run with different initial conditions, unknown equations of state and viscosities, 
untested breakup criteria, and varying chemistries, precludes robust or rigorous scien- 
tific conclusions being generated from comparison with data. But, despite the daunting 
modeling challenge that lies ahead, the analyses performed thus far have demonstrated 
a keen sensitivity between the bulk properties used in the calculations and experimental 
correlation analyses. 
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